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ABSTRACT 

Radiation interacts with matter via exchange of energy and momentum. When matter 
is moving with a relativistic velocity or when the background spacetime is strongly 
curved, rigorous relativistic treatment of hydrodynamics and radiative transfer is re- 
quired. Here, we derive fully general relativistic radiation hydrodynamic equations 
from a covariant tensor formalism. The equations can be applied to any three- 
dimensional problems and are rather starightforward to understand compared to the 
comoving frame-based equations. Current approach is applicable to any spacetime 
or coordinates, but in this work we specifically choose the Schwarzschild spacetime 
to show explicitly how the hydrodynamic and the radiation moment equations are 
derived. Some important aspects of relativistic radiation hydrodynamics and the dif- 
ficulty with the radiation moment formalism are discussed as well. 
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1 INTRODUCTION 



Radiation interacts with matter. The dynamics and energy 
contents of matter are influenced by surrounding radiation 
while radiation itself is concurrently determined by sur- 
rounding matter. So one needs to solve the hydrodynamic 
and the radiative transfer equations simultaneously and self- 
consistently. 

When matter is moving with a relativistic velocity, vari- 
ous relativistic effects should be taken into account: Doppler 
shift, aberration, time dilation, relativistic beaming, and so 
forth. These make physical quantities observer-dependent: 
one needs to specify carefully the observer from whose point 
of view physical quantities are measured. The same is true 
for matter and radiation in curved spacetime. Gravitational 
redshift, time dilation, loss cone effect, and so forth should 
be taken into account. To consistently deal with all these 
relativistic effects, one needs to build the radiation hydro- 
dyn amics within a f ully relativistic framework. 

iThomasl il93Clt) formulated special relativistic theory 
of radiative transfer to dea l with rad i ative viscosity in a 
differentially moving media. ILindauistl jl966l) subsequently 
developed a covariant form of general relativistic photon 
transport equation: a partial differential equation in time, 
space, angle, and frequency. A straightforward solution to 
such a fully angle- and frequency-dependent radiative trans- 
fer equation would be a direct integration of the multi- 
dimensional, partial differential equation. Although this 
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can be done in principle by a finite difference m ethod 
jHauschildt fc Wehrsell99ll : lLiebendorfer et ai]l2004h . char- 
acteristic methods which reduce the transfer equation to 
a single ordinary differential equation along characteristic 
rays are m ore practical a nd successfully applied to special 
relativistic iMihalasl 1980) and general relativistic problems 



relativistic (Mihalas 198U ) and general relativistic problems 
<Schmid-Burgldll97a [SchindeJll988t ISchinder fc Bludmanl 
19891: IZane et al.lll996t) . However, even these characteristic 
methods have been limited so far to spherically symmetric 
and steady cases. And if one is interested in the dynamics of 
the gas flow rather than in the details of the emergent spec- 
trum, simpler description of the radiation field may suffice. 

ILindauistl lll966tl also derived relativistic radiation mo- 
ment equations from the transport equation. In the ra- 
diation moment formalism, a finite number of angular 
moments, generated from the sp ecific in te nsityj are used 
to describe the radiation field. iThornel (Il98lf) general- 
ized the moment formalism to projected symmetric trace- 
free (PSTF) tensors to derive the relativistic version of 
moment equations up to an arbitrary order. PSTF for- 
malism has been successfully a pplied to a var ie ty of 
spherically symmetric problems ( Flammang 1982, 1984; 
Turolla fc Nobilil Il988t iNobili. Turolla. fc Zampieril Il99ll. 



1993tlZampieri. Turolla. fc TreveJll993t) . 

Almost all aforementioned formalisms are based on the 
comoving frame: matter and radiation quantities and their 
directional derivatives are defined and described within the 
comoving frame. It is most convenient to define various mat- 
ter and radiation quantities and describe their interactions 
in the comoving frame, a frame that moves along with mat- 
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ter. But, the directional derivatives in the comoving frame 
are rather awkward because the velocity of matter varies in 
space and in time. 

One can also present the radiation hydrodynamic equa- 
tions in Eulerian framework: derivatives appear in much sim- 
pler form in the Eulerian frame. The most straightforward 
Eulerian frame is a frame fixed with respect to the central 
object, which we call the fixed frame. The methodology of 
using comoving radiation quanti ties with i n fixe d-frame co- 
ordinates was first employed bv iMihalasl Jl980l) . which we 
will call the mixed-frame formalisiri TPeadTllr99^r also started 
from covariant tensor equations for energy-momentum to de- 
rive the mixed-frame radiation hydrodynamic equations for 
spherically symmetric systems: matter and radiation quan- 
tities and their derivatives were in covariant forms under 
fixed-frame coordinates while the interactions between radi- 
ation and matter were described by the comoving quantities. 
This approach makes the equations easier to interpret and 
apply to astrophysical problems. Recently, special relativis- 
tic radiation hydrodynamic equations for three -dimensiona l 
problems were derived by the same approach (|Parkll2004) . 
In this work, we extend this mixed-frame formalism to ar- 
bitrary spacetime or coordinates. We show specifically how 
the three-dimensional, general relativistic radiation hydro- 
dynamic equations are derived in the Schwarzschild geom- 
etry. We also discuss some aspects of relativistic radiation 
hydrodynamics as well as the difficulty in closing the radi- 
ation moment equations. We put c = 1 in most equations 
except in a few cases where c is explicitly written for clarity. 



2 TENSOR EQUATIONS 
2.1 Matter 

The energy-momentum tensor of an ideal gas is 



naff 



u g U a uP+P g g° 



(1) 



+ 



where U a is the four velocity of the gas and ui g ~- 
P g the gas enthalpy per unit proper volume which is the 
sum of the gas energy density e g and the gas pressure P g . 
The functional dependence of the gas enthalpy uj g on the 
gas temperature T can become comp lex, especially in the 
transrelativistic regime feervicdll986|) . 



2.3 Radiation Hydrodynamic Equations in 
Covariant Form 

Since mass and energy are equivalent in relativity, it is the 
particle number density, rather than the mass density, that 
is conserved in relativistic hydrodynamics. The continuity 
equation, therefore, should represent the conservation of the 
number density, 



(nU° 



= 0. 



(3) 



In the absence of any external force other than the ra- 
diative interaction, the total energy-momentum tensor of gas 
plus radiation is conserved, 

+ R oc0} 



IT 



= 0. 



(4) 



One can define the radiation four-force density to specif- 
ically describe t he energy and momentum transferred from 
radiation to gas iMihalas fc Mihalaslll984f) . 



G a = - I dv I dQ[ X I(n,u) -n]n a , 



(5) 



where x ls the opacity per unit proper length and r\ the 
emissivity per proper unit volume. The combinations r\jv 2 
and v\ are again frame-independent scalars. 

This radiation four-force density is equal to the diver- 
gence of the energy-momentum tensor for matter, 



nag 



G a , 



(6) 



and to the minus the divergence of the radiation stress ten- 
sor, 



R al 



-G a . 



(7) 



When micro-physical processes for the interactions between 
radiation and matter are known, equation J1J can be put 
into two separate equations © and @ and solved. 



3 SCHWARZSCHILD SPACETIME 

The above tensor equations can be explicitly writ- 

ten out in any spacetime or coordinates . In this work, we 
specifically choose the Schwarzschild spacetime to show ex- 
plicitly how to derive more familiar hydrodynamic and ra- 
diation moment equations and various transformation rela- 
tions among radiation moments. Much s imple r, spherically 
symmetric case has been shown in IParld lll993l) . 



2.2 Radiation 

The energy and momentum of the radiation field are repre- 
sented by the radiation stress tensor that consists of zeroth, 
first, and second moments in angle of the radiation field, 



7(n, v)n a nP dvdQ, 



(2) 



where I(x a ; n, v) is the specific intensity of photons moving 
in direction n on a unit sphere of projected tangent space 
with the frequency v measured by a fiducial observer, p a the 
four-momentum of photons, and n a = p a /hu. Since I„/V 3 
which is proportional to the photon distribution function 
and vdvdQ. are frame-independent scalars, R a/3 is a con- 
travariant tensor. 



3.1 Schwarzschild geometry 

The most familiar form of the Schwarzschild metric is 

dr 2 — —gapdx^dx 13 
dr 2 

= F 2 dt 2 T - r 2 (d0 2 + sin 2 8d(j> 2 ), 



(8) 



where M is the mass of the central object, m = GM/c 2 , and 

is the lapse function. 

The four-velocity of the gas, defined as 

dx a 



U 



dr 



(10) 
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satisfies the normalization condition 

U a U a = -1, (11) 
from which one can define the energy parameter 

y = -U t 

= [r 2 + ((7 r ) 2 +r 2 (rf/ e ) 2 +r 2 (rsin^) 2 ] 1/2 . (12) 



3.2 Tetrads 

It is not trivial to define physical quantities for matter or 
radiation in curved spacetime. However, one can always set 
up some tetrad, and since a tetrad frame is a locally inertial 
frame, physical quantities can be straightforwardly defined 
as in flat spacetime. Among various tetrads, fixed and co- 
moving tetrads are the most relevant ones. A fixed tetrad is 
an orthonormal tetrad fixed with respect to the coordinates 
and has base = d/dx 1 that can be expressed in terms of 
coordinate base d/dx' as 



dt 
d_ 
dr 
d_ 

d§ 
d_ 

dj> 



ld_ 

r dt 

dr 
ld_ 
rd9 
1 



(13) 



9 



r sin 6 8<j) 



Physical quantities defined in the fixed tetrad are those mea- 
sured by a fiducial observer who is fixed with respect to 
the coordinates. This observer sees matter moving with the 
proper velocity v, which is related to the four- velocity as 



v l = — = 



U a e- 



U* -U a et a ' 



(14) 



where U 1 and U* are the spatial and time parts of the fixed 
tetrad components of the four-velocity and e-" 1 is the a-th 
component of the tetrad base e;. Specifically, 



V 

y 

—rU , 

y 

-rAnOU^ 

y 



(15) 



where a and $ denote tetrad's bases. The components of 
Lorentz transformation A a s(v) are 



A f 3 



7 

JVj 

j + V V- 



(18) 



7-1 



We can also express the comoving tetrad in terms of 
coordinate base, which is useful when converting the tetrad 
components to and from covariant ones: 



9 
dLc 

9 
df cc 



7 9 , 3 19 19 

+ y Tv rTT + -yv e -— + ^ ; — 777, 

r at or r do r sin 6 d<p 



7 d „ 



_9 

dr 



d 



K< 1 v 2 r 89 y ' ' v 2 rsmOd^ 

r B dt + 17 ' v 2 dr 



(19) 



1 + (7"1)^ 



ld_ 

rd6 



+ (7-1) 

+ (7-l)^~ 
v 2 r 36 



^V r V<f, 9 

v 2 dr 



+ 



l + (7-l) 



r sin i 



Similarly, applying the inverse Lorentz transformation 
A a p(—v) yields the inverse transformation from the comov- 
ing tetrad to the coordinate base: 



1 

fdt 



die 



dr 



Otco 



1 + (7 - 1) 



- 7«^ 

9 
dfco 



9 
d<f>c, 



+ (7-!) — 5 1- (7- 1 ) — 5^^" 



Since v is a three vector defined in the tetrad, Vi = v l 
also define the Lorentz factor 7 for v as 



We 



7 = (1 -v 2 ) 



-1/2 



y 
r' 



(16) 



where v 2 = v ■ v — ViV 1 — v 2 + v 2 + v 2 . The value of v 2 at 
the horizon is always 1 regardless of IP: a fiducial observer 
fixed at the horizon always sees matter radially falling in 
with velocity c. 

A comoving tetrad moves with velocity v' relative to 
the fixed tetrad and therefore is related to the fixed tetrad 
by the Lorentz tranformation, 



9 a/3 r \ 9 
„ . = A &(v) — -, 



(17) 



ld_ 

rd9 



9 , i,v g v r 9 



l + (7-l)5 



+ (7 _ i)I!2H* _1 



9 



+ (7-1) 



« 2 90c 
- + (7- 

10 

v<t,ve 9 



7 "dLo +(1 ' v 2 9f c , 



(20) 
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3.3 Radiation Moments 

Radiation moments can be denned from the specific inten- 
sity, I v (x a ;n), either in a fixed tetrad or in a comoving 
tetrad as done in flat spacetime. 

The radiation energy density is equal to the zeroth mo- 
ment 



E 



I v dvdQ,, E co 



(21) 



where dQ and dQ co are the solid angle elements in the fixed 
and comoving tetrad, respectively. The radiation flux con- 
sists of three first moments 



IvtfdvdQ,, F l co = I Vco n co dv co dQ. co , (22) 



F l = 



where i — r,0, 4>. The radiation pressure tensor is symmetric 
and consists of six second moments 



Ivfin 3 dvdQ,, PH =11 I Vco n co n :i co dVcodQ. co .{22>) 



pi] _ 



From the definition of the radiation stress tensor (Eq. 
the tetrad components of the radiation stress tensor for 
the fixed tetrad are simply 



R 



&0 



( E 


F r 


F° 


p4> 


F r 


prr 


pre 


pr<j> 


F e 


prff 


prB 


pB<t> 


V Ft> 


pr<j> 


p8<t> 


p4"l> 



(24) 



For the spherically symmetric radiation field, F e = F^ = 
pre _ pr<t> _ q and p ee = p ^ = ( E _ p" )/ 2 . The comov- 
ing tetrad components are similarly 



R, 



a/3 





F' 

* CO 


F e 

' CO 


/.' 

' CO 


pr 

' CO 


prr 

CO 


prO 

CO 


prfi 

CO 


F° 

' CO 


prO 

CO 


pr9 

CO 


p84> 

± CO 


x CO 


prep 

1 CO 


p e<p 

1 CO 


T)4>4> 
1 CO 



(25) 



The contravariant components of the radiation stress 
tensor R a/3 are related to the fixed tetrad components by 
the transformation 



a Br a Br 13 - - 
p a P 

dx^ dx° 

and contain all the curvature and coordinate specifics, 



R*e = 



Y- 2 E 
F r 

p-l F e 



F r 

V 2pr 



Vr- 



El 



r sin 9 
r sin e 

p e<t> 

r 2 sin 9 

p4>4> 

r 2 sin 2 9 



(26) 



(27) 



n U r sin e r sin e 

Since tetrad components R a/3 and R^o are Lorentz ten- 
sors, they are related by the Lorentz boost A a p{v), 



R, 



A a ,(-v)A%(-v)R 



? A/5 



(28) 



Substituting equations 12411 and 1251 yields the stan- 
dard transformation law between the fixed-frame radia- 
tion moments and the comoving-frame radiation moments 
jMihalas fc Mihalaslll984l ; lMunier fe Weaverll98r3) . which is 
reproduced here for completeness: 



E co 
F i 



7 2 [E 
— , y 2 v l E + 7 



5} + ( 



(29) 



F J 



pi] 

CO 



2 i 7 m 

7 v v h — 7 



Ofc H 5— v v k 



(30) 



,7-1 i 



v 1 5t + v J 61 + 2 — = — v l v 3 Vk 
1^ 



F 



yZ yZ 



(31) 



Since this transformation is between tetrad quantities, it 
is valid in curved spacetime as well, as long as the proper 
velocity v is used. Replacing v with — v yields the inverse 
transf ormation. As discussed in detail in lMihalas fe Mihalasl 
dl984l) . the comoving- frame radiation flux can be signifi- 
cantly different from the fixed-frame flux due to the advec- 
tion of radiation energy density and pressure, v l E + VjP' 1 - 1 , 
and one should carefully determine which flux is the appro- 
priate one, especially in a lower-order approximation. 



3.4 Radiation Four-Force Density 

The energy and momentum transfer rates, measured by a 
comoving observer, are given by the comoving tetrad com- 
ponents of the radiation four-force density, 



G, 



dVco I dQ C o[XcoIv co - I)co]llco. 



(32) 



In terms of the usual heating function r co , the cooling func- 
tion A co (both per unit proper volume) and the mean opacity 
Xco (per unit proper length), all in the comoving frame, 



G C o — r co A C o 

(~<i ^ TP' 

° CO XcO-TcOl 



where 

Fco 
XcoF co 



I dv C o I d^lcoXco^Uco^co • 

c 



— I dVco j dQ.coXcoIv co , 

du co I dQcoVco, 



(33) 
(34) 

(35) 
(36) 
(37) 



Expression in the form of equation l|34|l is possible only if the 
thermally emitted (and/or scattered) photons are isotropic 
in the comoving frame so that the net momentum they re- 
move from the gas is zero. 

The components G a are related to the comoving tetrad 
components G™ by the transformation 



' 1 . / V 

/-iot OX ^jj 

— £~<~ r CO- 

dXco 

Substituting dx a /dXco from equation 11911 yields 



(38) 



G' 
G r 

G° 

G 4> 



G C o + ViGco 



G co + IVrGlo + ~ , ^ ' VrViGco 
V 

cto + JveGlo + 7 5 - veViGlo 



(39) 



1 



Gto + JV^Gto + 7 1 V<l,ViG co 
V 
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4 RADIATION HYDRODYNAMIC 
EQUATIONS 

4.1 Continuity Equation 

The continuity equation in the Schwarzschild cordinates JHj 



+ 



ld( sin enU 9 ) + -?-(nU' t -)=0, (40) 



which is the conservation equation of the particle number, 
rather than the mass density, because the mass density p in 
relativity necessarily contains the internal energy that is not 
always conserved. 

The continuity equation for spherically symmetric, 
steady-state flow becomes 



iirr nil 1 = — N = constant. 



(41) 



One should note that the velocity U r in equation I4H is the 
radial component of the four-velocity U a which is different 
from the proper velocity v T . It is the proper velocity v r = 
y~ 1 U r that is always equal to c at the black hole horizon 
(r = 2m) regardless of the value of U r while U r is not 
necessarily equal to c at the horizon. This applies to non- 
spherical cases as well. 



4.2 Hydrodynamic Equations 

The projection tensor PoP projects a tensor on a direction 
perpendicular to U a , 



Pj = gJ + U a U p = 8 J + U a U?. 



(42) 



Relativistic Euler equation is obtained by projecting equa- 
tion © with P a s to get Pp a T f3X ., x = Pp a G 13 , which becomes 

Lo g U a .,pU? + g af} P g , + U a U l3 P g ,p = G a + U a U G . (43) 

Radial part of Euler equation is obtained by fixing a = 



Trt dU r rrl dU r 
UJg U—+u Ja U — 



m 

i 

v 

r- 



+ ^[^(t/y-r-^ir) 2 ] 



— [(rU 8 ) 2 + (r sin 6»I7 ) 2 ] 

dt or ox 1 

= -yU r G* + [1 + r~ 2 (U r ) 2 ]G r + r 2 U r U e G B 

+r 2 s,m 2 6U r U' p G' l >. (44) 

The term mr~ 2 on the left-hand side represents the grav- 
itational acceleration and r 2 r~ 1 [(r(7 9 ) 2 + (r sin 9U' 1 ') 2 ] the 
centrifugal acceleration. The right-hand side of the equa- 
tion consists of all four components of G a , not just G r . In 
a sph erically symmetric case, the equation becomes JParkl 
Il993l) 

y dU r l d(U r ) 2 m yU r dP 9 y 2 dP g 

r dt 2 dr r 2 T Lo g dt Lo g dr 



V 



V 



X.coF co , 



(45) 



and the comoving- frame flux F£ a , rather than the fixed- 
frame flux F r , is directly responsible for the radial accel- 
eration. Therefore, even when the net flux F r measured by 
an observer at rest with respect to the coordinate is small, 
gas particles can still experience significant radiative force 
in the presence of significant radiation energy density E and 
pressure P 1J in the optically th ick medium, and slow down 
efficiently jPark fc Millerlll99ll) . 

The ^-part of the Euler equation is 

rrt dU B rn dU B 
dt " dx 1 



1 , 



<p\2 



+ 2cj g -U r U-ujgSmdcos8(U' p ) 

+ u B u^ + - ldP < 



dt 



dx 1 



yU u G l + V- 2 U r U B G r + [1 + T Z {U u f]G u 



+r 2 sin 2 e^U^G*, 



(46) 



and the </>-part 



idU* 



g „ _ + UgU _ 



(JgU 

+ 2Lj g -U r U 4 ' +2aj„coteU B U 4 ' 



4> T f 
4>r<t 



dt 



^Pg , rj<t> jji 9Pg 

' dd> dx* 



= -yU 9 G l +r-' 2 U r U 4, G r +r 2 U B U^G B 
+ [l + r 2 sin 2 f9([/ ) 2 ]G^. 



(47) 



Energy equation is more readily obtained by projecting 
Aion @ along the four- velocity U a T al3 . t/3 = U a G a , 

t 9 f U A „rr» 9 f^A , rrtdPg , , A dPg 



- nU t* (<*) _ nU *JL (<££) + u**i + ui^L 

dt V n J dx 1 \ n J dt dx 1 

= -yG l + F~ 2 U r G r + r 2 U B G B +r 2 sin 2 OU^G^ (48) 
_ a — v 

CO J1 CO * CO- 

Again, the exchange rate of energy between matter and radi- 
ation is expressed by heating and cooling functions measured 
in a comoving frame. 

4.3 Radiation Moment Equations 

The zeroth moment equation, i.e., the radiation energy equa- 
tion, is obtained by taking the time component (a = t) in 
equation J7J and expressing G a with the corresponding ra- 
diation moments, 



+ 



±d_E J_d_ 22 

r 2 dt F 2 r 2 dr K ' 

1 



d 



Yr sin 6 d6 y ' Fr sin 9 

= -G* = --^(r co -A co + XcoViF l co ) 



(49) 



In equation l|49|l . both the temporal and spatial derivatives 
are applied to the fixed-frame radiation moments while the 
interaction between matter and radiation is described by the 
comoving frame quantities. This 'mixed-frame' representa- 
tion simplifies the equation and makes each term easier to 
understand. For example, A co — r co is the energy input from 
the matter to the radiation field through the usual radia- 
tive heating and cooling processes and Xc ViFl is the work 
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done by the photons to the matter via transfer of mome ntum 
through scattering or absorption iPark fc Millerlll99lTl . 

When the spherically symmetric radiation field does not 
vary in time nor interacts with matter, equation 1491 be- 
comes the 'generalized luminosity conservation' equation, 



4nr 2 T 2 F r = L„ 



(50) 



where the constant Loo is the luminosity measured by an 
observer at infinity. A static observer at a given radius r 
finds the luminosity to be a function of radius, 



4TrrF r 



2m /r 



(51) 



This is the well-known gravitational redshift effect: locally 
measured (within the gravitational potential) fixed-frame lu- 
minosity can be larger than the luminosity seen by an ob- 
server far away (outside the gravitational potential) . For ex- 
ample, the fixed-frame luminosity near the compact object 
can be super- Eddington eve n when the luminosity measured 
far away is sub-Eddington l|Park][l992). 

The r-part of equation 10 yields the radiation momen- 
tum equation in r-direction, 



dF r 
dt 



, dP" 
dr 



r sin 9 89 



d , . opr es , r dP r 
sm9P H : — -— — 



+ ™ (E + P rr ) + — (2P rr - P 99 - P**) 
= -G r 



(52) 



= -YxcoF r co - r 7 u r (r co - a co ) - r^—^-vrv.xcoF^. 

In relativistic treatment, w 1 -order term v r (T co — A co ) and 
the redshift term (m/r)(E + P)/r should be included as well 
as the correction factors F and y. These relativistic terms 
are negligible in slowly moving, optically thin flow. How- 
ever, when tv > 1, they significantly affect t he radiation 
field as well as the dynamics of the gas flow ||Par klll99Ct 
iPark fc Mineilll99lt |Parklll992l : iPark fc Milleilll995h . 

The radiation momentum equation in ^-direction from 
equation Q is 



r dt r dr v I rsinfl 89 



d I ■ aT -,ee\ 
1 sin 9P 



1 dP 9 * 2T pre 



1 



r sin 9 dcj> 
= -rG e 



r tan t 



(53) 



-XcoF co - 7«e(r co - A co ) 



and that in rA-direction is 



IdF* Id ( ^ ldP°* 
r dt r dr K > r 89 

?Hp r <P , 2 pO<t> 
r r tan 9 

-rsin6'G" A 



-veViXcoFlo, 



dd> 



+ 



(54) 



-V^ViXcoFlo. 



-XcoFf - 7^(r co - A co ) 



In a static non-relativistic case, equations Q52H54J re- 
duce to the familiar diffusion equation 



F l = - -r— P tJ . 

Xco 



(55) 



4.4 Discussion 

Although reducing the angle-dependent radiative transfer 
equation to the radiation moment equations lowers the di- 
mensionality of the problem from seven (space, time, angle, 
and frequency) to four (space and time), it comes with a 
price: while the number of moment equations (14911 and 152t 
I54H is only four, the number of all radiation moments to be 
determined, E, F l , and P tJ , is ten. This is a generic problem 
with moment equations: equations at each order always con- 
tain successively higher order moments. Astronomers rou- 
tinely employ some kind of closure relation to close the equa- 
tions, for example, the Eddington factors. In complex radi- 
ation hydrodynamic calculations, an educated guess of the 
Edddington factors a s functions of the optical depth is so me- 
times tried (see, e.g.. lTamazawa et aljfl975l lParklll990fl . 

But it can't be stressed enough that the correct form 
of the Eddington factors or any other closure relations can 
only be obtained by solving the f ull angle-dependent ra- 
diative transfer equation (see e.g., lAuer fc Mihalas! Il97fj. 
iHummer fc Rvbickil Il97ll lYin fc Milled Il995l) . which has 
been rarely tried fo r three-dimensional flows and much 
less for relativistic (ISchinder fc Bludmanl [l989) or time- 
dependent flows . Even in one-dimensional systems, the vi- 
able method, iterative variable Eddington method, requires 
expensive computations and, moreover, when the flow mo- 
tion is not monotonic various angle-frequency components 
couple along a ra y and even posing boundary conditions be- 
come a problem iMihalas!ll980 | : [Nobili. Turolla. fc Zampieril 
Il993f) . Iteration method bv ISchmid-Burgkl il978i) seems 
to be the only workable approach for the general rel- 
ativistic system, but it works only for spherically sym- 
metric cases. Although there had been a suggestion for 
the general shape of the Eddginton factors for arbi- 
trary three-dim ensional radiation hydrodynamic systems 
jMinerbolll978ft . its validity has yet to be verified by true 
angle-depend ent radiative transfer calcu lations for a vari- 
ety of cases dSchinder fc Bludmanl H989T) . Needless to say, 
in special cases where the specific intensity or the radiation 
moments can be directly calculated, the incomplete radi- 
ation moment equati ons can be eschewed altogether (e.g., 
IChattopadhvavll2005h . 

Other related problem with the radiation moment for- 
malism is that the exact nature of the radiation field can be 
properly described only by the infinite number of moments. 
Using a finite number of moments by terminating higher or- 
der moments or assuming some kind of closure relation can 
be a good enough approximation in certain conditions, but 
not necessarily so in any situation: it can lead to spurious 
pathological behaviours such as a ra diation shock, especially 
when the velocit y gradient is steep llTurolla fc NobiiHll988l : 
|Pullemondlll999l) . 



5 SUMMARY 

General relativistic, radiation hydrodynamic equations are 
derived using a covariant tensor formulation. Matter and 
radiation quantities are defined in the fixed and comoving 
tetrads and transformed to corresponding covariant forms. 
The interactions between matter and radiation are described 
by the radiation four-force density defined in the comov- 
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ing tetrad and transformed to the covariant form. Con- 
servation of energy and momentum separately for matter 
and radiation yield the hydrodynamic and radiation mo- 
ment equations. Current approach is applicable to any three- 
dimensional systems in any spacetime or coordinates, but in 
this work, we show explicitly how the equations are derived 
for the Schwarzschild spacetime. Certain aspects of relativis- 
tic radiation hydrodynamics and the fundamental limitation 
of the radiation moment formalism have been discussed as 
well. 
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